Implications of the Reduction Principle for 
Cosmological Natural Selection 

Lee Altenberg 

Associate Editor of BioSystems; Ronin Institute; altenber@hawaii.edu 

Smolin (1992) proposed that the fine-tuning problem for parameters of the Standard Model 
might be accounted for by a Darwinian process of universe reproduction — Cosmological Natural 
Selection (CNS) — in which black holes give rise to offspring universes with slightly altered 
parameters. The laws for variation and inheritance of the parameters are also subject to CNS if 
variation in transmission laws occurs. This is the strategy introduced by Nei (1967) to understand 
genetic transmission, through the evolutionary theory of modifier genes, whose methodology is 
adopted here. When mechanisms of variation themselves vary, they are subject to Feldman's (1972) 
evolutionary Reduction Principle that selection favors greater faithfulness of replication. A theorem 
of Karlin (1982) allows one to generalize this principle beyond biological genetics to the unknown 
inheritance laws that would operate in CNS. The reduction principle for CNS is illustrated with a 
general multitype branching process model of universe creation containing competing inheritance 
laws. The most faithful inheritance law dominates the ensemble of universes. The Reduction 
Principle thus provides a mechanism to account for high fidelity of inheritance between universes. 
Moreover, it reveals that natural selection in the presence of variation in inheritance mechanisms 
has two distinct objects: maximization of both fitness and faithful inheritance. Tradeoffs between 
fitness and faithfulness open the possibility that evolved fundamental parameters are compromises, 
and not local optima to maximize universe production, in which case their local non-optimality 
may point to their involvement in the universe inheritance mechanisms. 
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I. SMOLIN'S COSMOLOGICAL NATURAL 
SELECTION HYPOTHESIS 

Smolin [1] has proposed that the free parameters of 
the Standard Model — and of the laws of nature, more 
generally — may be the result of a process in which grav- 
itational singularities cause new universes to form with 
slightly different laws of nature, and the entire branch- 
ing ensemble of universes comes to be dominated by laws 
which produce the greatest number of gravitational sin- 
gularities, through a process of 'cosmological natural se- 
lection'. The elements of the hypothesis are (slightly 
paraphrased from [2]): 

1. The world consists of an ensemble of universes. 

2. Black hole singularities bounce and evolve to initial 
states of expanding universes. 

3. Hence there is a fitness function equal to the aver- 
age number of black holes produced by a universe 
initiated in such a bounce transition. 

4. At each such creation event there is a small change 
in universe properties leading to a small random 
change to the average number of black holes pro- 
duced. 

5. Under very mild assumptions for the fitness func- 
tion, the ensemble converges after many universe 
creation steps to one in which almost every mem- 
ber is near a local maximum in terms of universe 
production. 



6. Therefore, for a randomly chosen member of the 
ensemble, almost every small change in the parame- 
ters of the Standard Model from its present value ei- 
ther leads the production of black holes unchanged 
or leads to a decrease. Since our universe can be 
assumed to be randomly chosen we conclude that 
if the hypotheses made above are true, almost no 
change in the parameters of the Standard Model 
from the present values will increase the numbers 
of black holes produced. 

I do not address here the plausibility of any of these 
elements with respect to the physics. Critiques and re- 
sponses in this regard are discussed in [2-11]. I also do 
not address the unsettled issue of how probabilities over 
multiverse spaces should be interpreted in light of the da- 
tum of our universe. Rather, I address the evolutionary 
dynamics. 

In recent years, the theory of evolutionary dynamics 
has broadened considerably beyond the phenomena orig- 
inally identified by Darwin — adaptation and diversifica- 
tion. Augmentations of the 20th century "Modern Syn- 
thesis" of Darwinian natural selection with Mendelian ge- 
netics have grown into a body of theory extensive enough 
to be called an "extended synthesis" [12]. Among the ear- 
liest developments in this area was theory to account for 
the evolution of the genetic mechanisms themselves. 

Smolin's proposal leaves mostly open the matter of uni- 
verse 'genetics' — the mechanisms that might produce 
changes to the Standard Model parameters from universe 
creation events. Some discussion is presented in [10]: 
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The hypothesis that the parameters p change 
by small random amounts should be ulti- 
mately grounded in fundamental physics. We 
note that this is compatible with string the- 
ory, in the sense that a great many string 
vacua likely populate the space of low-energy 
parameters. It is plausible that when a re- 
gion of the universe is squeezed to the Planck 
density and heated to the Planck tempera- 
ture, phase transitions may occur, leading to 
jumps from one string vacua to another. But 
so far there have been no detailed studies of 
these processes which would have checked the 
hypothesis that the change in each generation 
is small. One study of a bouncing cosmol- 
ogy in quantum gravity also lends support to 
the hypothesis that the parameters change in 
each bounce [13]. 

One need not wait for these issues to be resolved before 
investigating how universe inheritance laws would enter 
into the evolutionary dynamics. Moreover, one can ask: 
What if universe inheritance laws are themselves subject 
to evolution? Here, another point made by Smolin [14] 
is relevant: 

[T]hc evolution of laws implies a breakdown 
of the distinction between law and state. An- 
other way to say this is that there is an en- 
larged notion of state — a meta-state which 
codes information needed to specify both an 
effective law and an effective state, that the 
effective law acts on. 

The 'meta-state' to account for the evolution of universe 
inheritance laws would be information coding for the in- 
heritance laws, which allows them to vary between uni- 
verses and thus be subject to CNS. 

II. MODIFIER THEORY FOR THE 
EVOLUTION OF INHERITANCE MECHANISMS 

Evolutionary theory in biology took just this approach 
in trying to understand the evolution of biological inheri- 
tance mechanisms when Nei [15] posited the first modifier 
gene model. Mendel's laws of genetic transmission had 
hitherto been handled as exogenous parameters in the 
mathematical models of evolution of the Modern Syn- 
thesis. There was no reason to believe, however, that the 
material basis for Mendel's laws was any different from 
the material basis for any other character of organisms. 
The development of molecular biology confirmed that the 
mechanisms of inheritance — recombination, mutation, 
DNA repair, etc. — utilized the same objects as organ- 
ismal traits involved with survival and fitness: proteins, 
nucleotides, regulatory sequences, gene interaction net- 
works, and self-organizing structures and activities in the 
cell and organism. 

Nei's modifier gene concept [15] explicitly "enlarged 
notion of state" with "a meta-state which codes infor- 



mation needed to specify both an effective law and an 
effective state." The meta-state took the form of a gene 
posited to control the recombination rate between two 
other genes under natural selection. 

The first rigorous stability analysis of this model by 
Fcldman [16] showed that new forms of the modifier gene 
could grow in number if and only if they reduced the 
recombination rate when they occurred in populations 
near a stationary state. The same Reduction Principle 
[17] was observed for modifiers of mutation rates and 
dispersal rates as well as recombination (see [18-24] for 
a sample of the early studies). 

The Reduction Principle proposes that populations at 
a balance between transmission and natural selection will 
evolve to reduce the rate of error in transmission. The 
fundamental implication is that evolution has two differ- 
ent properties that it operates to maximize: organismal 
fitness, and the faithfulness of transmission. 

The mathematical basis of the Reduction Principle was 
discovered by Karlin [25] (although he did not realize 
it) in a fundamental theorem on the interaction between 
growth and mixing: in a system of objects that (1) are 
changed from one state to another by some transforma- 
tion processes, and (2) grow or decay in number depend- 
ing on their state, then greater mixing produces slower 
growth or faster decay. To be precise, Karlin's theorem 
states: 

Theorem 1 (Karlin [25, Theorem 5.2]). Let M be an 

irreducible stochastic matrix, and D a diagonal matrix 
with positive diagonal elements. If D ^ c I for any c e 
R, then the spectral radius r({(l—a)I+aM)D) is strictly 
decreasing in a G [0, 1]. 

This phenomenon is not restricted to finite state or 
discrete time models. Karlin's theorem has now been ex- 
tended to general resolvent positive operators on Banach 
spaces [26] . This relationship between mixing and growth 
is therefore an extremely general mathematical property. 

The implication for the theory of cosmological natural 
selection is that the parameters of the Standard Model 
may be optimized not simply for universe replication 
rate, but also for faithfulness of universe replication, and 
these two objectives may not coincide, but could poten- 
tially involve tradeoffs. Therefore, any finding that pa- 
rameters of the Standard Model are not optimal for the 
rate of universe replication, via the production of black 
holes in the case of Smolin's proposed mechanism, cannot 
be taken on face value as falsification of the Cosmological 
Natural Selection hypothesis. A further examination of 
how such parameters may be constrained by the mecha- 
nism of universe inheritance, or how they may impact the 
faithfulness of universe replication, must be considered. 

Main Implication. Parameters in the Standard Model 
that are not local optima for the rate of production of 
black holes may be involved in the inheritance mecha- 
nisms of universe creation and subject to contravening 
selection for the preservation of universe properties. 
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III. ILLUSTRATION WITH A BRANCHING 
PROCESS MODEL 

The modifier gene methodology may be applied to in- 
vestigate the evolution of hypothesized universe inheri- 
tance mechanisms through addition of another degree of 
freedom to the universe state to allow for variation in the 
inheritance laws. 

Variation in the transmission laws can be drawn from 
a very large choice of alternatives. Variation in the mag- 
nitude of changes in the free parameters of the Standard 
Model is a natural choice. However, to model the con- 
sequences of variation in magnitude one has to include a 
model of the map from parameters to reproductive prop- 
erties, a speculative and complicating exercise. The other 
principle alternative is to vary the probabilities among a 
fixed set of changes in state. The simplest variation in 
probabilities is to scale equally all of the transitions be- 
tween states. Such variation is employed here. 

The process of universe creation through black hole 
creation is naturally modeled as a continuous time 
branching process with an infinite number of types and 
with aging to account for individual universe develop- 
ment. The purpose here, however, is merely to illustrate 
how universe inheritance mechanisms can evolve, and for 
this it is sufficient to analyze a simplified branching pro- 
cess with discrete time, finite number of types, and no 
aging. The multitype Galton- Watson process is defined 
as follows. 

First, some notation is introduced. The whole num- 
bers are represented as No := {0, 1,2,.. .}. Vectors are 
represented as x,v,e, etc. and matrices as N,P,D,I, 
etc.. Matrix elements are ATy = [N]ij. The vector of 
ones is e, its transpose e T , and is the vector with 1 at 
position % and elsewhere. Let D x := diag[xi] represent 
a diagonal matrix with diagonal elements Xi. A column 
stochastic matrix P satisfies P > and e T P = e T . An 
irreducible nonnegative matrix A has for every pair 
some t > 1 such that [A*]ij > 0, and if it is aperiodic, 
A t > for some t>l. Letting \i(A) be the eigenvalues 
of A, the spectral radius is r(A) — maxi |Aj(A)|, and the 
spectral bound is s(A) = maxj Re(A,(A)), where Re(x) 
is the real part of x. 

Let x £ Nq be the n-vector of the number of indi- 
viduals of each of n types. The probability that a sin- 
gle individual of type j produces a vector of offspring 
numbers x is p(x,j). Starting with an initial popula- 
tion x(0) at time t = 0, each individual produces an 
offspring vector x according to p(x,j), and these vectors 
are summed to create the population at the next time 
step, x(l). The reproduction process then repeats. This 
defines the branching process {x{t) : t = 0, 1, 2, . . .}. 

The expected number of offspring of type i from parent 
3 is 

N ij = x iP( x J)- 

It is assumed that exists for alH, j — {1, . . . , n}. The 
expectation for the process is E[£c(t)|a;(0)] = iV*a;(0). 



The following theorem gives the asymptotic behavior 
of a super-critical branching process. 

Theorem 2 (Kcstcn and Stigum [27] from [28, Theorem 
1, p. 192]). 

Let x(t) £ Nq, t £ No be vector of the number of indi- 
viduals Xi of type i in an n-type Galton- Watson branching 
process. Assume the following: 

1. The process is nonsingular, i.e. it is not true that 
for every j there is an i such that p(et,j) = 1. 

2. The process is supercritical, i.e. r(N) > 1. 

3. N is irreducible and aperiodic. 

Let v(N) be the Perron vector of N (the eigenvector 
associated with the spectral bound of N) normalized so 
e T v(N) = 1. 

Then (almost surely, a.s.) 

lim = cv(N), 

t^oo r(N)* 

where c > is a random variable such that Pr[c > 0] > 
if and only if for each i,j £ {1, . . . , n} 

L V : = Y X i l °g( X i)P( X >j) < 00 • I 1 ) 

Let us now frame a multitype branching process for 
cosmological natural selection. Let U represent the set of 
possible universe types, with a finite number \U\ of possi- 
ble types. Time is censused in epochs t = 0, 1, From 

one epoch to the next, each universe produces offspring 
universes according to the multitype Galton- Watson pro- 
cess. The persistence of a universe between time epochs 
is accounted for by considering it one of its own offspring. 

Some structure is now introduced for the expectation 
matrix N . We decompose in terms of the total ex- 
pected number of offspring and the distribution among 
that total. Let 

be the expected total number of offspring from type j. 
This value j3j represents the fitness of universe type j. 

Let Tij > be the expected fraction type i among the 
offspring of j, so for each j, J2 ieU Tij = 1. Then N is 
decomposed into a fitness component and an offspring 
distribution component: = Tij/3j. Define the matrix 
T := [Tij]ij e n. Then in matrix notation: 

N = TD P . 

The offspring distribution T is further decomposed. The 
expected proportion of offspring that differ from parent j 
is given as fij . Thus 1 — fij is the expected fraction of type 
j's offspring that are also type j, i.e. Tjj = 1 — fij. The 
expected fraction of offspring of type j that are type i 
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among the non-type-j offspring is Pij, hence T^- = HjP%j 
for i j. This decomposition in matrix form is 

T = (I D fx )+PD li = I+(P (2) 

where T and P are stochastic matrices. It should be 
noted that many different probability measures p(x,j) 
can produce the same matrix of expected proportions T. 

Now let us include variation in the branching pro- 
cesses. We enlarge notion of state i G U to a meta- 
state (a, i) G U' := T x W, where T is the set of pos- 
sible transmission laws for U. A transmission law for 
the branching process is the probability measure p a {x,j), 

p a : Nq W ' xU —> [0, 1], that defines the branching process. 
An important simplifying assumption made here is that 
the transmission law states a G T are themselves trans- 
mitted faithfully in offspring universes. A model in which 
the transmission laws themselves transform in universe 
creation is necessary to a more complete theory, but is 
deferred to further study. By making the transmission 
law faithfully replicated, we know that its evolution is 
due solely to its effects on the universe replication be- 
havior. 

Further, variation in transmission law states a G T is 
assumed for now to be neutral, in that it leaves unchanged 
the vector of the mean number of offspring, (3. 

Variation in inheritance laws will be encapsulated 
solely through a parameter m a G [0, 1] (a G T) that 
scales equally the expected proportion of offspring uni- 
verses that differ from their parents' type. In the popu- 
lation genetics literature this is referred to as linear vari- 
ation [29] . In matrix form, the scaling by m a enters into 
the expected proportion matrix T„, a e T, from (2) as 

T a = (I- m a D^) + maPDf, =I + m a (P- I)D„. 

Thus the matrix of the expected numbers of offspring 
for universes within the transmission class a G T is 

N a = [I + m a (P — I)Dfj]Dp. (3) 

Remark 1. A model in which there is variation in the 
magnitude of changes between parent and offspring uni- 
verses would also be manifest as changes in N, but they 
would not be linear variation. For example, let changes 
in a universe parameter jj follow a Gaussian distribu- 
tion (a case considered by Vilenkin [11]), parameterized 
by m a , to give expected distribution of offspring i having 
parameter -fi, 

e -(li -ljf/ m l 

In linear variation, m a scales all (m a ) equally for i ^ j, 
and disappears from the derivative of T t j (m a ), but clearly 
this is not the case in (4). The analysis of r(T(m a )Dp) 
under such non-linear variation is very much an open 
problem. 



Now let us examine the evolution of the multiverse 
ensemble 

x := (xi,x 2 ,...,x v ) G Nq X|W| 

where there are v different transmission laws in the en- 
semble, where each transmission law a G T is carried by 
sub-population x a G Ng W '. The matrix of expected num- 
bers for the entire ensemble is a direct sum of diagonal 
blocks N a : 

N := N 1 ®N 2 ®-- -®N V . 

The diagonal block structure is the result of the trans- 
mission laws themselves being faithfully transmitted. 
With framework defined, we have the following result. 

Theorem 3 (Domination of the Most Conservative 
Transmission Law). 

1. Consider an ensemble of universes reproducing ac- 
cording to a branching process, where transmission 
law a G T determines the transmission law parame- 
ter m a G (0, 1] which scales the expected proportion 
of universes i G U that differ from their parent uni- 
verse types j <ElA. 

2. Let the matrix of expected numbers of offspring uni- 
verses of type i from a universe of type j be 

N a = [I + m a (P - IjD^Dp, for each aeT, 

where fi is the vector of the fractions of offspring 
that are not identical to their parents, m a scales fi, 
[3 > is the vector of expected numbers of offspring 
universes for each type in U, and P is the matrix of 
expected distributions of offspring type i from parent 
type j ± i, Pa = 0. 

3. In addition make the technical assumption that for 
each a G T using p a (x, j) in (1) that Lij < oo for 
all i,j G U. 

4- Assume that the ensemble branching process is 
supercritical, so the spectral radius is r(N) = 
max aer r(N a ) > 1. 

5. Let stochastic matrix P be irreducible and aperi- 
odic, and let some f3i ^ Let the transmission 
parameters be ordered so that mi <m 2 < ■■ -ra v . 

Then 

= (ci«(iVi),0,...,0) a.s., 

where c\ is as in Theorem 2. 

Therefore, the ensemble of universes becomes dom- 
inated by the transmission law with the highest ex- 
pected preservation of universe types in universe creation 
events. 
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Proof. For the branching process for each transmission 
type a G T, Theorem 2 gives 



lim 



T x a (t) = c a v(N a ) a.s.. 



i^oo r(AT a )* 
Suppose that r(7V b ) < r(N a ). Then 
1 . r(N b y 



l 



= c 6 t,(JV 6 ) lim 3^=0 a.s.. 
t-i-oc r{I\ a ) z 

Thus we need to know when r(N b ) < r(N a ). This occurs 
when rnh > m a , as shown by Theorem 1. An equivalent 
form from [26] is that 



dm 



s(D + mA) < s(A), 



(5) 



where A is irreducible and has nonnegative off-diagonal 
elements, and s() refers to the spectral bound, which is 
the spectral radius if the matrix is nonnegative. For the 
case here, A = (P — ^D^Dp, so r(A) = since e T (P — 
^D^Dp = ( e T — e T )D fJ ,Dp = 0, showing e T to the be 
left Perron vector of A. Thus r([I + m a (P - IjD^Dp) 
is strictly decreasing in m a . 

Therefore, if m a < m b then r{N a ) > r(N b ). Hence 
mi < m,2 < . . . < m v gives r(N) = r(Ni) > r(N 2 ) > 
... > r(N v ). Thus 



t^co r(N) t 



for 6 e {2,3,..., v}. 



This produces the claimed result that 
1 



-(aji(i),a;2(i), • ■ -,x u {t)) 



t^,oo r(Ny 

= (ci»(JVi),0,...,0) a.s.. 



□ 



Remark 2. It should be noted that no simplifying as- 
sumptions are made regarding the relationship between 
the expected number of progeny universes [3j and the 
universe inheritance distributions . No constraints are 
made on the (3j other than they be positive, and that T 
be irreducible. The result here is completely general in 
this regard. 

In this branching model, we see that the transmission 
law with the highest preservation of the laws of physics 
during universe creation becomes the only typical mem- 
ber of the ensemble as time goes on. 



reproduction rates, we may give type a its own vector of 
fitnesses (3 a , and then the matrix of expected offspring 
numbers becomes 



N a = [I + m a (P — I)Dfj]Dp a 



(6) 



It is clear that if m a < nib then there is some region of 
(3 a values that maintains r(N a ) > r(AT h ), and this range 
includes values for which (3 ai < f3 bi for each universe type 
i G U. 

Hence, universe type (a, i) can become dominant over 
another type (b, i) even though (3 ai < /?(,$, due to m a < 
m b . 

Therefore, to find a universe type (b,i) G T x U with 
fitness greater than our own universe type (a, i) does not 
necessarily falsify Cosmological Natural Selection, but 
rather points to the involvement of property a in uni- 
verse transmission and suggests that transmission law a 
is more conservative than transmission law b. 



IV. ADDITIONAL LIMITS TO OPTIMIZATION 

Another obvious caveat for predictions from Cosmolog- 
ical Natural Selection is the assumption that all of the 
parameters of the Standard Model can be varied in any 
direction. Since the origin the parameters is mysterious, 
and the hypothesis that they vary as a result of universe 
reproduction is speculative, the possibility that the pa- 
rameters are constrained to vary only in some subspace 
can't be ruled out, in which case optimality arguments 
have to be made with respect to variation within that 
subspace and not for a single variable. Varying a single 
variable might appear to increase universe reproduction, 
but not be allowed within the space of variation. The 
problem of constraints in variation has a long history in 
evolutionary theory [30] , as well as the general question of 
whether evolution optimizes organismal phenotypes [31]. 

Even in the case where universe creation varies the 
fundamental parameters in all directions, an additional 
key point is raised by Vilenkin [11]: evolution produces 
a distribution over the parameter space when it reaches 
a mutation-selection balance, and the intensity of muta- 
tion pressure compared to selection pressure determines 
the closeness of a typical universe to a fitness optimum. 
Hence, falsifiable predictions cannot be made without 
some knowledge of the mutation process. The question of 
how concentrated or spread out the distribution is around 
optimal types under mutation-selection balance is a prin- 
cipal issue in the subject of quasispecies theory [32, 33]. 



III.l. When Transmission Fidelity is in Conflict 
with Universe Production 

The variation in universe transmission laws is assumed 
to be neutral as framed Theorem 3: ra a varies without 
varying the expected offspring production ft. However, 
nature need not produce such a clean separation of ef- 
fects. If the transmission type a G T also affects the 



V. DISCUSSION 

Theorem 3 shows that the ensemble of branching uni- 
verses comes to be dominated by universe transmission 
laws which are the most faithful in preserving the fun- 
damental constants among offspring universes. This is a 
manifestation of the Reduction Principle from population 
biology applied here in a new setting. 



6 



The ubiquitous manifestation of this phenomenon is 
due to the fundamental mathematics of how mixing in- 
teracts with heterogeneous growth. The most general 
characterization of this interaction is for operators on 
Banach spaces in [26] , that 

4-s{mA + V) < a(A), 
dm 

where A is resolvent positive operator on a Banach space, 
V is a real-valued operator of multiplication, sQ is the 
spectral bound, and m > 0. Resolvent positive operators 
A include Schrodinger operators, second order elliptic op- 
erators, diffusions, integral operators, and in general, any 
generator of a positive strongly continuous semigroup. 
This theorem derives from Kato's [1982] generalization 
to Banach spaces of Cohen's [1981] theorem that for fi- 
nite matrices A, s(A + D) is convex in diagonal matrices 
D. 

This dynamic has two principle implications for the 
Cosmological Natural Selection hypothesis: 

1. It gives a causal mechanism to incorporate universe 
transmission laws into the framework of cosmolog- 
ical natural selection, and predicts that transmis- 
sion of the laws of physics between parent and off- 
spring universes should be very conservative. 

2. It demonstrates that in addition to selection for the 
maximal number of offspring produced, selection is 
also indirectly minimizing the probability of change 
between parent and offspring. Should there be a 



tradeoff between these two properties, then the ob- 
servation of non-optimal values of Standard Model 
parameters with respect to black hole production 
may point to the involvement of those parameters 
in the universe transmission law. 

The model analyzed here, while more general than a 
toy model, also leaves out of consideration how trans- 
mission laws may be involved in their own transmission. 
There are also many other higher-order phenomena in 
evolutionary dynamics that may prove valuable to in- 
clude in models of cosmological natural selection, includ- 
ing the evolution of mutational robustness, the evolution 
of evolvability, and the evolution of modularity in the 
genotype-phenotype map. 

While the Cosmological Natural Selection hypothe- 
sis remains largely speculative, inclusion of higher-order 
phenomena of evolutionary dynamics may provide ad- 
ditional guidance in deriving falsifiable predictions from 
the theory. 
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